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Introduction générale 


Les phénomènes vibratoires jouent un rôle déterminant dans presque toutes les 
branches de la physique: mécanique, électricité, optique, acoustique, etc. Malgré 
leur grande diversité, ils sont régis, en tout cas dans le domaine linéaire, par les 
mêmes lois de comportement et peuvent être étudiés au moyen du même outil 
mathématique. L'homme s'est intéressé aux phénomènes vibratoires lorsqu'il a 
construit les instruments de musiques. 

Malgré les connaissances acquises par les anciens, il faut attendre le début du 
17e siècle pour que Galilée (1564-1642) démontre que le ton d'un son est 
déterminé par la fréquence des vibrations. 

Très longtemps, on a étudié les vibrations des machines et des structures presque 
uniquement dans le but de les atténuer et, si possible, de les supprimer. Cette 
préoccupation est encore essentielle mais n’est pas la seule. On construit 
actuellement de plus en plus de machines ou d’appareils qui utilisent les 
vibrations mécaniques pour remplir la fonction désirée. Les vibrations sont 
parfois perturbatrices et doivent être combattues dans plusieurs domaines par 
exemple : 

- Les machines ou certains organes de machine sont une cause 
d’imprécision, de bruit, d’usure prématurée et de fissure, entraînant 
finalement la rupture de la pièce. 

- Les vibrations des voitures, des avions, des trains ou des bateaux 
provoquent, en plus des inconvénients précédents, l’inconfort des 
voyageurs et diminuent parfois la sécurité de conduite de ces véhicules. 

- Les vibrations des grandes structures métalliques peuvent prendre, dans 
certains cas, des proportions catastrophiques. 


Dans ces chapitres, nous allons étudier le comportement des systèmes avec un ou 
plusieurs degrés de liberté. Nous allons nous limiter à des systèmes où les 
équations de mouvement sont des équations différentielles linéaires. Ce qui nous 
permet de décrire de diverses caractéristiques importantes de vibrations. Vous 
apprendrez à analyser les vibrations libres et forcées avec ou sans 
amortissement. 
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CHAPITRE I 


Oscillations libres non amorties : Système à un degré de liberté 


1.1 Généralités sur les vibrations 

1.1.1 Mouvement périodique : 

Définition : C’est un mouvement qui se répète à intervalles de temps réguliers, cet intervalle est appelé 
période (T) qui s’exprime en seconde (s). 

Pour les mouvements rapides, on utilise la fréquence : / exprimée en Hertz (HZ) 

1.1.2 Mouvement vibratoire : 

Définition : Un mouvement vibratoire est un mouvement périodique se produisant de part est 
d’autre d’une position d’équilibre. On peut aussi définir un mouvement vibratoire par sa 
fréquence/. La fréquence indique le nombre d’oscillations complètes (dans le sens aller retour) 
se produisant par seconde. 

On peut établir la relation entre la fréquence et la période : 



La période T des oscillations est le temps mis par le système pour revenir à une 

position identique quelque soit le choix de cette position. C’est aussi, le temps mis 

pour faire une oscillation complète ou un « aller-retour ». 

Mathématiquement, le mouvement périodique de période T est défini par: 

A tout instant t, x(t + T) = x(t) 

1.1.3 Mouvement vibratoire libre 

Définition : les vibrations libres sont les vibrations qui résultent lorsqu’on écarte un système de 
d’équilibre ou on lui donne une vitesse initiale, puis on le laisse vibrer librement. 

Exemples : Une masse accrochée à un ressort - un pendule simple - le balancier d’une horloge - 
d’un moteur tournant à vitesse constante etc. 

1.1.4 Mouvement vibratoire sinusoïdal 

Définition : un mouvement vibratoire est sinusoïdal, si un point vibrant possède une élongation du type : 

y(t) — A sin (mt -I- <p) 

> La grandeur y(t) est appelée l’élongation (ou la position) 
à l’instant t, l’élongation maximale ou l’amplitude du 
mouvement, elle varie entre -A et +A. 

> La quantité O) est la pulsation du mouvement et exprimée 
en ( ra< ^ / s ). 

> La quantité (o)t + <p) est la phase instantanée, exprimée en 
(radian, sans dimension), 

> l’angle (p est la phase initiale, correspond à la phase à l’instant t = 0. 

1.2 Vibration harmonique 

Définition : On appelle vibration harmonique tout système dont le paramètre x(t) qui la caractérise est une 
fonction sinusoïdale du temps : x(t) — A cos (&)t + cp) 


x(t) 



t(s) 


sa position 
la rotation 
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> La fonction cosinus est une fonction périodique de période 2n. Si T est la période temporelle du 
mouvement, on aura donc : 

[co(t + T) + cp] — [eût + cp] — 2n => ù)T — 2 n 

2 TC 

On en déduit l’expression de T en fonction de la pulsation : T = — 

Ù) 

> La fréquence/, nombre d’oscillations par seconde correspond à l’inverse de la période T:f = VT. 

Il existe d’autres expressions équivalentes pour la fonction x(t). En effet, la fonction sinus est équivalente à 
la fonction cosinus décalée de td/ 2. On peut donc écrire : 

x — A cos (tut + cp) — A sin (mt + <p) 

Donc : 

Les grandeurs caractéristiques d’une vibration harmonique sont : 

L’amplitude A, 

2n 

La période T, to — — — 2nf ; co: pulsation ,/: fréquence. 

La phase (p. 

1.2.1 Coordonnées généralisées d’un système physique 

Définition : Les coordonnées Généralisées sont l’ensemble de variables réelles indépendantes ou liées 
permettant de décrire et configurer tous les éléments du système à tout instant t. 

Par exemples : 

> un point matériel libre dans l’espace peut être déterminé par 3 coordonnées généralisées (x, y, z); 

> un corps solide peut être déterminé par 6 coord, génér. : 

• 03 coordonnées relatives au centre de gravité; 

• 03 coordonnées liées aux angles d’Euler (cp, xp, 8). 

^ Les coordonnées généralisées d’un système de P points 

matériels et Q corps solides sont défini par : N — 3P + 6Q 

On note : 

Les coordonnées généralisées : qi(t),q 2 (t), 9w(0- 

Les vitesses généralisées : qi(t), q 2 (t), q' w (t). 

1.2.2 Degré de liberté 

Définition : Le degré de liberté est le nombre de coordonnées généralisées indépendantes, nécessaires pour 
configurer tous les éléments du système à tout instant : d — N 

Où, le nombre de coordonnées généralisées liées, pour configurer tous les éléments du système à tout instant 
moins (-) le nombre de relations reliant ces coordonnées entre elles : d — N — r 
d: Degré de liberté ; 

N : Nombre de coordonnées généralisées 
r : Nombre de relations reliant ces coordonnées entre elles. 

Exemples : 

■ Un disque de masse m et de rayon r, roule sans glisser sur un plan horizontal. 

Ici on a deux coordonnées généralisées x et 6 donc N — 2 . 
x et 8 sont liées avec une relation: x — r8 donc :r — 1 . 

Le nombre de degrés de liberté d — N — r — 1. — 



/. 
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■ Un système mécanique constitué de 02 points matériels et M 2 reliés d’une tige de longueur 1. 
M 1 (x 1 ,y 1 ,z 1 ) : 3 
M 2 (x 2 ,y 2 ,z 2 ) : 3 


N — 6 


1 


L’équation de liaison : l — ^ (x x — x 2 ) 2 + (y 1 — y 2 ) 2 + {z 1 — z 2 ) 2 = c te Mi 

=> r — 1 => d = 5 


M 2 


1.3 Equation différentielle du mouvement 

Dans ce cours, on établi l’équation différentielle en utilisant le formalisme de Lagrange. L’intégration de 
cette dernière permet de donner l’équation du mouvement. 


1.3.1 Formalisme de Lagrange 


Ce formalisme repose sur la fonction de Lagrange (L — T — U ). L’ensemble d’équations du mouvement 
s’écrit : 



• L : Fonction de Lagrange ou Lagrangien 

• T : L’énergie cinétique du système; 

• U : L’énergie potentielle du système ; 

• q t : est la coordonnée généralisée et q t est la vitesse généralisée du système. 
Pour un système à un degré de liberté, (N= 1 ou ddl=l) l’équation du mouvement s’écrit : 


d /5L\ / dL 

dtvdô/ \dq 


= 0 


Remarques : 


• Pour un mouvement unidimensionnel x, l’équation de Lagrange s’écrit : 


d /ôi\ /ÔL\ 

dt VaxJ VdxJ 


= 0 


• Pour un mouvement rotationnel 9, l’équation de Lagrange s’écrit : 


d 

dt 



1.3.2 Exemples d’oscillateurs harmoniques 


Exemple 1 : Pendule élastique vertical 

Un pendule élastique est constitué d’une masse suspendue à un ressort de 
raideur k et peut donc osciller verticalement avec une élongation x(t). 

Le système nécessite une seule coordonnée généralisée x(t) qui peut décrire 
le mouvement de la masse m et de l’extrémité mobile du ressort. 

Donc le système a un seul degré de liberté d=N=l. 
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• L’énergie cinétique du système: T — -mx 

• L’énergie potentielle du système: l’énergie U emmagasinée dans le ressort dépend de 
l’allongement des 2 extrémités du ressort. Elle s’exprime: 

U = \k(x 2 - Xi) 2 = \kx 2 avec x 2 = x ; x 1 = 0 


X 

dU — F r ■ dx = — (— kxdx ) => U — j kxdx — -/ex 2 


La fonction de Lagrange : L — T — U 
L'équation de Lagrange : £ g) - (g) = 0 


L = imx 2 — -kx 2 
2 2 


dL 

= mx : 

dx 

dL , 

— - —kx 

dx 


d ( dL \ - 
dt \dxJ 


mx 


mx — (— kx) — 0 


On divisant par m x -\ — x — 0 

r m 

k / k 

Le rapport — étant positif et en posant : ù) 0 — — on obtient l’équation différentielle d’une vibration 


harmonique de la forme : x + co 0 2 x — 0. 




La pulsation co 0 ne dépend que de la masse m et de la raideur k du ressort, est appelée « la pulsation 


propre » du système. 

•S La masse oscille donc indéfiniment avec une période propre T 0 donnée par la relation suivante: 



Exemple 2 : Pendule pesant simple 


Un pendule simple est constitué d’un solide de petite dimension de masse m suspendu à un point fixe O par 
un fil inextensible de longueur L. Ecarté de sa position d’équilibre, il oscille dans le champ de pesanteur 
terrestre g. 


Les coordonnées du système : 

( x - 
ml 
ty = 


= l sin 6 =» x — W cos 6 
y — l cos 6 => ÿ = —W sin 6 


■ L’énergie cinétique du système : T — - mv m 

T — ^m(x + ÿ) 2 = im/ 2 0 2 (cos d 2 + sin d 2 ) 

=> T = -ml 2 è 2 
2 
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■ L’énergie potentielle du système : U — mgh (h est la hauteur de m par rapport à un plan de 

référence donnée.) 

NB : On a deux possibilités pour calculer la valeur du déplacement h, selon le choit de l’origine des 
énergies potentielles (U(0)=0), ce choit doit avoir lieu lorsque la masse est dans sa position d’équilibre 0 = 0. 
L’énergie potentielle correspond à l’énergie potentielle de pesanteur. 


1 er cas : si on choisi comme origine des énergies potentielles l’axe ( Ox ) on a donc : 
h — — l. cos 9 (Le signe moins vient du fait que la masse m est inférieur à l’axe choisi). 
Dans ce cas : U — —mgl. cos 0. 

2 eme cas : si on choisit comme origine des énergies potentielles (U(0)=0) l’axe (O'x 1 ). 

À l’équilibre, on aura : h =1 — 1. cos 9. Dans ce cas : 


U - mgl(l — cos 0) 

Calcule du lagrangien : L = T — U 

1 er cas : On remplaçant T et U dans L on trouve : L = -ml 2 9 2 + mgl. cos 9 


L’équation de Lagrange : 


d /5L\ / dL 

dt \dé) \d9 


f d (dL\ .. 

— r = m. I. 2 9 (L 

dtVd0/ ' 

i (dL\ 

(jg) = —m.g. I.sin9 ( 2 ) 

v. (1) — (2):m. Z. 2 9 + m.g.l.sin9 = 0 


( sin 9 ~ 9 

Dans le cas des faibles oscillations, les angles sont très petits on a : j cqs g ~ ^ — 

On aura donc m. I. 2 9 + m.g. 1.9 = 0, En divisant par m. I. 2 on trouve : 


1 


0 +j0 = 0 

C’est l’équation de l’oscillateur harmonique de pulsation propre : ù) 0 = 

On trouve enfin : 0+m o 2 0 = 0. 



2 eme cas : L = ^ ml 2 9 2 + mgl( 1 — cos 0) 


d /5L\ 

dtvàé/ 



= m. I. 2 9 


—m.g. l.sin 9 


■( 1 ) 

( 2 ) 


(1) — (2): m. I. 2 9 + m.g. l.sin 9 = 0 
On aura donc m. I. 2 9 + m.g. 1.9 = 0 , On divisant par m. I 2 on trouve : 

6 + j 6 = 0 Avec w 0 = j , et on retrouve bien le même résultat. 


Université Ferhat Abbas-Sétif- 


F acuité de technologie 


Tronc commun ST 


N. AKLOUCHE 



CHAPITRE I Oscillations libres non amorties Système à un degré de liberté 

1.3.3 Solution de l’équation différentielle du mouvement 

L’équation différentielle (EDF) du mouvement est de la forme : 

x + ù) 0 2 X = 0 

C’est une équation différentielle du second ordre sans second membre dont la solution sous la forme 
complexe est de la forme : x(t) = Ae at 

La dérivée première de la fonction x(t)) (la vitesse) : x(t) — Aae at . 

La dérivée seconde de la fonction x(t) (l’accélération) : x(t) = Aa 2 e at 
On remplace dans l’EDF : Aa 2 e at + œ 0 2 Aae at — 0 => Ae at (a 2 + n> 0 2 ) = 0 

Or Ae at ^ 0 => a 2 + co 0 2 = 0 donc a — ±jco 0 
Donc la solution aura la forme: x(t) = A 1 e JÙJot + A 2 e~ JÙJot 
Selon la relation d’Euler : e ±JÙ)()t — cos <n 0 t ± jsin co 0 t 

-> x(t) = A 1 (cos co 0 t + jsin œ Q t) + A 2 (cos m 0 t — jsin co 0 t) 
x(t) = ( < A 1 + A 2 ) cos m 0 t +J(A 1 — A 2 )sin co 0 t = C cos co 0 t + Dsin co 0 t 
Tel que : C =(A 1 + A 2 ) et D = j(A 1 — A 2 ) 

Donc x(t) = C cos co 0 1 + Dsin co 0 t est aussi une solution de l’équation différentielle. 

Si on pose : C — a. cos d et D — a.sin 9 , on aura : x(t) = a. cos 0 cos co 0 t + a.sin 6 sin co 0 t 
cos(x — y) = cos x cos y + sin x sin y donc : x(t) = a cos( m 0 t — 9) — a cos( m 0 t + <p) , cp — 9 + ^ 

donc : x(t) = a cos( m 0 t + <p)Tel que : a — sic 2 + D 2 et 0 — arc tang(j) 


1.4 La force dans le mouvement harmonique 

1.4.1 Exemple du pendule élastique vertical 

C’est le cas d’une masse m accrochée à l’extrémité libre 
d’un ressort et se déplaçant sans frottement suivant une 
direction Ox vertical (voir figure). 



Ressort à vide Equilibre avec Mouvement 
une masse 

(a) (b) (c) 
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A l’équilibre : il y a deux forces qui agissent sur la masse m ; son poids et la force de rappel du ressort 
tension due au ressort : 

£ F = P + T = O* => mg - kàl = 0 

• P : Poids de la masse m. 

• T : Force de rappel du ressort. 


En mouvement : La deuxième loi de Newton (principe fondamental de la dynamique), nous permet 
d’écrire : 

F — my 


Pour un système à une dimension : 


d 2 x 

F — m — mx 
cLt 2 


Après projection on obtient mx — mg — k(x + Al). En utilisant la condition d’équilibre précédente on 
obtient : 

mx = —kx => mx + kx — 0 


Or : co 2 n = — => mx + mœ 2 n x = 0 

u m u 

=> mx — —kx — —mœ 2 0 x ; C’est la force de rappel due au ressort 

avec k - ma ) 2 0 = cte 

Donc, la force dans les mouvements harmoniques simples est proportionnelle et opposée au déplacement 
et constitue une force de rappel. 


1.4.2 L’étude d’une vibration harmonique en termes d’énergies 

Nous voulons montrer que l’énergie totale (mécanique), E=T+U, est constante et déduire la valeur de cette 
constante. Pour cela prenons x = Acos(œ 0 t + <p), alors : 


111 1 
E — T + U — -mx 2 + — kx 2 - -mA 2 co 2 0 sin 2 (ù) 0 t + <p) + -kA 2 cos 2 (co 0 t + cp) 


Sachant que (sin 2 tp + cos cp 2 — 1) et en utilisant la relation : k — mœ 2 0 

alors, E = -mA 2 a) 2 n — - kA 2 — constante 

7 9 u 9 


T IJ 

1 max u max 

Nous retrouvons ici le fait que l’énergie mécanique de ce système ne varie pas. L’énergie totale est 
constante. 

On a: T — -mx 2 — -mA 2 ù) 2 0 sin 2 (co 0 t + cp) — -mA 2 ù) 2 0 [ 1 — cos 2 (m 0 t + cp)]. 


T =imA 2 m 2 0 [l-^] = ^ma) 2 0 [A 2 - x 2 ]. 


l x = 0 => T = T max =^mœ 2 0 A 2 . 


Si 


x = => T min = 0 
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D’autre part : U — -kx 2 — -ma) z 0 x 


fx — 0 


Si 


U — U min — 0 (position d'équilibre) 


x — +A 


U Umax 2 k* 0 


La figure suivante montre la variation des énergies cinétique, potentielle et totale en fonction de x : 



■S L’énergie se transforme d’une énergie cinétique à une énergie potentielle. 

•S Quand l’énergie cinétique diminue l’énergie potentielle augmente et vis versa. Cette propriété de est 
appelée conservation de l’énergie totale du système. 


1.5 Systèmes équivalents 

Définition : C’est un système simple qu’on représente en générale par un ressort équivalent ou une masse 
équivalente. 

1.5.1 Masse équivalente : Cas d’un ressort de masse non négligeable. 

m : La masse du ressort. 

Au repos : 

• Z : La longueur du ressort. 

• dm : masse élémentaire située à une 
distance y du point de suspension. 

En mouvement : 

• x(t) : Déplacement instantané de l’extrémité mobile du ressort. 

• dy : Déplacement de la masse élémentaire = y x(t) => sa vitesse = yx(t) 

o La masse linéique du ressort à une distance Z:p=y=>m = pZ. 
o La masse de l’élément dy du ressort : m s = p dy = y dy 

=> L’énergie cinétique = L toutes les énergies de ses éléments; 



^T = { 1 -(jdy).^m 2 ') 
(0 
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1 m 
2 y 3 


m 


T = ~ (~ )x(t) 2 <=> -m eq x(t) 2 m eq = — 


1.5.2 Ressorts équivalents : On a 3cas : 

1 er cas : Ressorts en parallèles (en oppositions) 


//////////// 


777777 



Ou 



////// 


<=> k 



J *(0 


777777 

L’élongation de chaque ressort est égale à x(t) donc : M. g — (/q + k 2 )x — keq.x 


k eq = fri + k 2 


2 eme cas : Ressorts en séries : 

Soit x t : l’élongation du ressort /q tel que : M . g — /q x t 
Soit x 2 : l’élongation du ressort k 2 tel que : M. g — k 2 x 2 

1 1 

=> x = x i + x 2 = M . g (— + — ) 

Ki K 2 


77/7/y 


l _ l l 


k± k 2 


'■eq 

léme 


/q /q 


et? fci + k 2 



777777 


< [> ^eq 


Xi 


x 2 



3 cas : Barre liée à 02 ressorts (Distance non négligeable ) 

L 

◄ ► 



<^> 



k en = 


(a+by 


q fc2 

fe 2 kl 


Si a=b, on aura : k eq = Jq + k 2 
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1.6 Analogie entre le système mécanique " Masse-ressort" et le système électrique "L-C". 


Système mécanique 

Système électrique 

Déplacement : x(t) 

Charge électrique q(t) 

Vitesse : x(t) 

Courant électrique i — 

Accélération : x 

Variation du courant : q 

Masse : m 

Inductance, bobine, self : L 

Ressort k 

Inverse de la capacité 1/C 

Force de rappel : k x 

d.d.p entre les bornes d’in condensateur : ^ 

Force d’inertie : mx 

d.d.p entre les bornes de la bobine : L q 

Energie potentielle : - k x 2 

Energie électrique : ~Q 2 

Energie cinétique : ^ mx 2 

Energie magnétique : ~kq 2 


CCI Points clefs 

Oscillations libres non amorties 
1. Pendule élastique vertical (m, k, x): 


x + o) 0 2 x = 0 - 

s 

2. Pendule pesant simple (m, l, G) 


( 

' 

k 

x(t ) = A cos (o ) 0 1 + tp) avec 

“ 0= ) 

— ( pulsation propre) 

V 


2n [in 

r 0 = — = 2 it - 

co 0 si k J 



f 

ù) 0 = ( pulsation propre) 

9 + ù) q 2 9 = 0 

0(0 = A cos ( o) 0 t + <p) avec 

< 

2n 

T 0 = =2 n 

[ "o J 

lï 

9 


> L’équation de Lagrange pour un mouvement unidimensionnel x ; ^ (jf?) ~ — 0; 

> L’équation de Lagrange pour un mouvement rotationnel 6 : ^ = 0; 

> L’énergie mécanique se conserve : T + U = Constante; 

171 

> Masse équivalente (Cas où la masse m du ressort n’est pas négligeable) : m eq — 

> Ressorts équivalents : 


f Ressorts k t , k 2 , ... k n en parallèles: k eq = k t + k 2 + 

111 

Ressorts k t , k 2 , ... k n en série : - — = — + — + 


, +k r 


+ r- 


Barre liée à 02 ressorts ( Distance non négligeable ): k eq = 


( a+b ) 2 

a^ + b^ 
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CHAPITRE II 

Oscillations libres amorties : Systèmes à un degré de liberté 

Introduction : Le pendule élastique comme le pendule pesant, se comporte comme un oscillateur 
harmonique à la condition de négliger tout frottement. Il oscille alors théoriquement sans jamais s’arrêter. 
En réalité, la masse se déplace dans un fluide (en général l’air) où il existe toujours des forces de 
frottement de type visqueux. L’oscillateur est alors amorti et fini par s’arrêter. 

II. 1 Oscillations libres amorties 

La présence de frottements implique une dissipation d’énergie sous forme de chaleur ; on observe alors 

• soit des oscillations dont l’amplitude diminue au cours du temps, 

• soit un retour à l’équilibre sans oscillation. 

On parle alors d ’ amortissement . L'expression de la force de frottement visqueux est la suivante : 

F q = -aq 


Tel que : 

a : est le coefficient de frottement visqueux, a : [Ai. s/m], 
q : la cordonnée généralisée du système ; 
q : La vitesse généralisée du système. 


Le signe moins (-) vient du fait que cette force s'oppose au mouvement en agissant dans la direction et le 
sens contraire à la vitesse. 

Dans un mouvement unidimensionnel x la force s’écrit sous la forme : 

/ = - a v = - axu 

II.2 Equation de Lagrange dans un système amorti 

En tenant compte de la force de type frottement fluide (coefficient de frottement visqueux a ), l’équation 
de Lagrange dans ce cas devient : 

d édL\ dL 

dt \dq ) dq Q 

Sous l’action des forces de frottements, le système dissipe (perde) de l’énergie mécanique sous forme de 
chaleur, il ya donc une relation entre la force F q et la fonction de dissipation D d’un côté et la fonction 
de dissipation et le coefficient de frottement visqueux a : 



et D — — aq 2 


L’équation de Lagrange dans le cas d’un système amorti devient : 


d 

dt 



dL 

dq 


dD 

dq 


II.2.1 Equation différentielle : Système masse-ressort-amortisseur 

Reprenons le cas du pendule élastique (vertical par exemple). L’étude de l’oscillateur amorti se fait de la 
même façon que précédemment mais en ajoutant la force de frottement visqueux. 

dD 


A une dimension, l’équation de Lagrange s’écrit : ^ = 


dx 


L’énergie cinétique du système : c’est l’énergie cinétique de la masse m : T= - mx 2 


L’énergie potentielle du système : c’est l’énergie emmagasinée dans le ressort U= - kx 2 
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La fonction de dissipation : 
La fonction de Lagrange : L 


= T-U 


D — - ax 2 
2 

L — - mx 2 — - kx 2 
2 2 


d fdL\ 

— —= mx 

dt \dxj 

dL , 

— = — kx 

dx 

dû 

— — ax 

V dx 


En remplaçant dans l’équation de Lagrange on aura : 


mx + kx — —ax 


.. a . k _ 

->> x + — x + —x = 0 


m m 


//////////// 


□ a 


m 


kx ax 

▲ à 


m 

▲ 


mx 


•S C’est l’équation différentielle du mouvement dans le cas d’un système libre amorti. 

•S Par rapport aux oscillations libres non amorties, on reconnaît un nouveau terme (£■ x) provenant 
de la dissipation d’énergie. 

S La forme générale : q + — q + —q — 0 

° m m 

•S Souvent l'équation différentielle est écrite sous une forme dite réduite : q + 2S q + U) 2 0 q — 0 


5 = — [l/Sl: Facteur d'amortissement. 

2 m L ' J 


Tels que : \ 5 

' f = — (Sans unité) : Rapport d’ amortissement. 

t “o 

À une dimension la forme réduite s’écrit : x + 28 x + &) 2 0 x — 0 

II.2.2 La solution de l’équation différentielle : Système masse-ressort-amortisseur 

L’équation différentielle du mouvement : x + 28 x + co 2 0 x — 0 

Il s’agit d’une équation différentielle linéaire du second ordre à coefficients constants sans second 
membre. 

La fonction x(t) = De rt est une solution particulière de cette équation différentielle à condition que r 
soit une des deux racines r 1 et r 2 de l’équation du second degré, appelée équation caractéristique. 

r 2 + 2 ôr + coq = 0 

La solution générale de l’équation prend la forme : x(t) = C x e r > 1 + C 2 e r ' 2t 
(r ± = -S + J S 2 - o ) 0 2 

Tel que: 1 ; On voit bien que la solution dépend des valeurs de 5 et w 0 . 


= -8-4^ 


(On 


l ei cas : ô < Mo (0 < <( < 1) : système sous- amorti ou faiblement amorti 
=> Ô 2 — ù) 0 2 < 0 


-1 

r i = S + Jj 2 (oj 0 2 

[r 2 = -ô-Jj 2 (œ q 2 

"a = V“0 

2 - ô 2 = w 0 V 1 - Ç 2 

T - 271 - . 

2n 

2n 

U (O a V w 0 2 -ô 2 

s 2 


Ü)Q 

1 2 

"0 


T 0 


/w 2 


Donc : T„= 


t o 


a 


; T a : pseudo-période 
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La solution : x(t) — C e 5t sinf<n n t + <p) 

Remarques : 

• x(t) représente un mouvement vibratoire. 

• L’amplitude C e~ St est décroissante : x(t) 
tend vers 0 quand t augmente. 

• l’élongation x(t ) va osciller en restant comprise 
entre - C e~ St et C e~ St .Ces deux exponentielles 
représentent l’enveloppe du mouvement de 
l’oscillateur c’est-à-dire les positions extrémales 

prises par x lorsque le temps s’écoule. 



I 


2 éme cas : S — oio = 1 ) : Amortissement critique : ry = r 2 = —ô 

La solution : x(t) = ( C 1 + C 2 t ) e~ St 

Si 5 — — et co 0 = — => a — a c — 2 s/ km : Valeur 

2m u v|m c 

critique du coefficient de frottement. 

Remarques : 

• x(t) n’ est pas oscillatoire car il ne contient pas un 
terme sinusoidal. 

• x(t) tend vers 0 sans oscillation quand le temps 
augmente. 

• Le système revient à sa position d’équilibre le plus 
rapidement possible. 

3 éme cas : 8 > cûo (f > 1) : système sur- amorti ou fortement amorti 
fr x = -S + Jô 2 - co 0 2 
{ r 2 = -8 - yjô 2 - co 0 2 



V> 


La solution : 

x(t) = e- St (C 1 e V« 2 -^o 2 t + c 2 e-#W‘) 

Remarques : 

• x(t) tend vers 0 sans oscillation quand le temps 
augmente. 

• x(t ) est un mouvement non sinusoïdal 



II.3 L’oscillateur harmonique électrique 

Nous allons voir maintenant qu'il existe un autre type d'oscillateur harmonique amorti dans un autre 
domaine de la physique : l'électricité. 

Soit un circuit électrique, constitué des 3 éléments de base misent en série : 

• un résistor de résistance R ; 

• un condensateur de capacité C ; 

• et une bobine d'inductance L. 
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Selon la loi de de Kirchoff : 
u R + u c +u L — 0 => Ri(t ) + ^ q + L ^ — 0 

= fl9+;<î + i? = 0 

/? l 

g + — g + - — g = 0 => g + 25 q + 6o 2 0 c7 = 0 

L LC 


1 

[ 5 = — 

1 

{ 8 = — 
2 L 

avec ( 

, 2L i Donc : 

q + 28 q + ù) 0 q = 0 => ( 

[T 

1 

£ 

O 

h 

II 

SI' 

1 



Remarque : 

• Pour un amortissement critique 8 — co 0 




Donc : 



II.4 Décrément logarithmique 

Définition : C’est le logarithme du rapport de 2 amplitudes successives des oscillations amorties. 


D 


= in 


*(£i) 

x(t 2 ) 


t 2 - ti + T a 


Où x(t 1 ) et x(t 1 + T a ) représentent les amplitudes des oscillations aux instants t { et (tj -1- T a ): 
généralement ces deux instants sont choisis comme correspondant à deux extrema successifs de même 
signe. Cette quantité mesure la décroissance des amplitudes pendant une période. 


x(L) x(L) 

D = ln— = ln . , 

x(t 2 ) x(ti + T a ) 

Pour un système amorti : 

x(t) = C e~ St sin(m a t + tp) 

^ = ln C e~ st ^ smjùig^ + (p) 

^ n C e~5(ti+r a ) sin^Cti + T a ) + (p) 

D = ln(e 5r “) = ÔT a 



ST a = Sj£==fr 0 


7^t t=2 ’ T ™ :donc: 

D = l n «h} = STa = 2 K 7 J= 


X(t 2) 


Remarques : 

• Pour plusieurs périodes : T — nT a ; t 2 — + nT a 




n i X (tl) i 

D — ln — — - = ln 

xÇt 2 ) 


xfti) nf 

7 " , = n8T a = 271—== 

Xity + nTa) ^/l - , 


La pseudo-période et le décrément logarithmique n’ont de sens que si le régime est 
pseudopériodique. 


II.5 Facteur de qualité (Facteur de surtension) 


Pour décrire l'amortissement d'un système oscillant mécanique ou électrique on emploie le facteur de 
qualité Q définit par l’expression suivante : 
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Q = 2ir^2f 

V |AE| 

• Emax '■ est l’énergie maximale stockée dans le système. 

• | AF | : est l’énergie perdue par cycle. 

• la notion de ‘qualité’ pour caractériser l’oscillateur, comme la grandeur qui traduit l'aptitude du 

système considéré à garder son énergie tout en oscillant. La qualité est d’autant meilleure que le 
£■ 

rapport est grand. 

II.5.1 Calcule du facteur de qualité : système masse-ressort-amortisseur (m. k. al 

Prenons l’exemple d’un système masse-ressort-amortisseur ( m , k, a ) faiblement amorti dont la solution 
de l’équation différentielle est sous la forme : 

x(t) = x 0 sin(£o 0 t + tp), x 0 — C e~ ôt et to a - Jco 0 2 — 8 2 « m 0 


On a d’une part : E max =^mœ 2 0 x 0 2 (Cf. chapitre I). 

D’autre part: A E — J^ +ta F(t)dx 

Tel que : F(t): est la force de frottement visqueux : F(t) = —a x(t) 


t+T n 


t+T n 


t+T n 


=> A E — j —a x(t)dx = J —a x{t)[xdt\ — —a j x 2 

t t t 

On a : x(t) = x 0 sin(o> 0 t + q)) => x(t) = x 0 m 0 cos (o> 0 t + <p) 


dt 


t+T n 


=> A E — —ax 2 0 œ 0 2 j cos 2 (m 0 t + cp) dt 

rt+T a r 

I cos 2 (œ 0 t + cp)dt — I 


“ 1 + cos 2(m 0 t + <p) 1 

dt => A E — — -aT a œo Xq 


2n 


T a — — ; ù) a ~ (jôq -+ AF — — ancooXQ 


(J0n 


o On retrouve bien une variation négative de l’énergie c’est-à-dire une perte d’énergie au cours du 
temps. 

o L’énergie perdue se transforme en énergie thermique ou elle se disperse en se diffusant dans le 
milieu avoisinant. 

En remplaçant dans l’expression de Q, on trouve : 

_ 2 n ° X ° _ tnw 0 Q __ mcüp _ <o 0 _ 1 

TtaXQÙ} 0 a a 28 2% 

II.5.2 Calcul du facteur de qualité : système électrique (RLC ) 

Dans un système mécanique (m, k ): E max =^mco 2 0 x 0 2 (cf. Chapitre I) 

Dans un système électrique (RLC): E max — ^ Lu) 2 0 q q 

Dans un système mécanique (m, k, a): AF = — naxQœ 0 => Dans un système électrique (RLC): 
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A E — —uRq^oiQ 

Donc : 


Remarque : 

• Plus l’amortissement est faible, plus la qualité du système oscillant est grande. Or Q est d’autant 
plus grand pour un œ 0 donné, que l’amortissement est faible. Un système très amorti a un Q 
faible. 



Q = 2n 


2 Lù) 2 0 ql 
nRqlo ) 0 



CQ Points clefs 
Oscillations libres amorties 

❖ L’équation de Lagrange pour un mouvement unidimensionnel x : 

d dL dD 1 -, 

— (— ) - — D = - a x 2 

dt K dx J dx dx 2 

(s = £- t 

♦> L’équation du mouvement : x + 2 S x + U ) 2 0 x = 0 =» < s 

( ^ w 0 


♦♦♦ x(t) = Ae rt est une solution particulière tel que : 


r x = —S + yjô 2 - vv 0 2 
r 2 = -ô-yjô 2 -w 0 2 


Amortissement faible ■ S < co 0 (0 < ^ < 1) S 2 — a) 0 2 < 0 
> La solution : x(t) — C e~ St sin(m a t + (p) 

co a = V^o 2 —S 2 : Pulsation des oscillations amorties. 
T a = — : pseudo — période. 

CO a 


avec 


Amortissement critique : ô = (Oo (£ = 1 ) 

> La solution : x(t) — ( C 1 + C 2 ) e~ St avec a — a c — 2 s/km 

Amortissement fort : S > to 0 > 1) 

> La solution : x(t) — e~ St (D t e^ s2 ~ w ° 2t + D 2 e~^ s2 ~ u> ° 2t ) 


î 

2Ç 


❖ Le Décrément logarithmique : D = ÔT a 

l'énergie stockée 


❖ Facteur de qualité : Q — 2 n.— 
1 I é 


ma) q 


Loi o 


énergie perdue par cycle 
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CHAPITRE III 

Oscillations forcées amorties : Systèmes à un degré de liberté 


Introduction : On a vu que l’amortissement des oscillations était dû à une diminution de l’énergie 
mécanique sous forme de chaleur dissipée. Pour compenser ces pertes d’énergies et entretenir 
(conserver) les oscillations, il faut une source d’énergie à travers une force extérieure. On va donc 
rajouter une force extérieure souvent dite excitatrice. 

Il va donc y avoir une force supplémentaire qu’il vaut mieux qu'elle soit colinéaire au mouvement et 
qu'elle soit le plus possible dans le sens du mouvement. 

Dans ce chapitre, on étudie la réponse d’un système amorti à 1 ddl à une excitation harmonique 
sinusoïdale produite par une force extérieure au système. Ce type d’excitation se rencontre fréquemment 
dans l’industrie (machines tournantes, ventilateurs, moteurs, pompes ...). 

III. 1 Equation différentielle du mouvement 


„ , , • . . d ,dL. dL dD „ 

a) Pour un mouvement de translatron, on ecrrt : — (— ) = + r ext 

dt dcj dcj dc[ 


o 


d f dL -, 
dt ^dè J 


dL 

d0 


dD 

dé 


Tel que : M {F ext ) = F ext . L = — \F ext \ 

M(F ext ) : Est le moment de la force appliquée [N.m]. 

Le moment : caractérise la capacité d’une force à faire tourner un objet autour d’un point. 

L : Le bras du levier : est la distance droite d’action de la force, 
r : La distance parcourue par la masse dans la direction de l’action de la force. 

III.l.l Exemple : système masse-ressort-amortisseur W/////////////A 

Reprenons le cas du pendule élastique (vertical par exemple). 

L’étude de l’oscillateur amorti se fait de la même façon que 
précédemment mais en ajoutant une force extérieure (F ext ). 

A une dimension, l’équation de Lagrange s’écrit : 

d ,dL. dL dD 

dt y dic' dx dx ext 

Prenons une force sinusoïdale appliquée à la masse m : F ext — F 0 sin cet. 

L’énergie cinétique du système : c’est l’énergie cinétique de la masse m : T = ^ 

L’énergie potentielle du système : c’est l’énergie emmagasinée dans le ressort : U = - kx" 

La fonction de dissipation 



ext 


t=0 



kx ax 


m 


n 


DIX 


ext 


D — -ax 2 

2 


2 --kx 2 
2 


9- 


mx 


La fonction de Lagrange: L =T-U= - mx 

rd / dL 

dt Vdx 
dL 

— = -kx 
dx 

dD 

— — = ax 
v dx 

En remplaçant dans l’équation de Lagrange on aura : mx + kx — —ax + F 0 sin cet 
On divise alors par m et on trouve : x + — x + — x — — sin ù)t 

m m m 


F 0. 


Souvent l'équation différentielle est écrite sous la forme réduite : x + 28 x + ù) 2 0 x — ^sin cet 
ô = —[1/51: Facteur d'amortissement. 

-Tl J 2m L ' J 

Tels que : < § 

' ^ = — (Sans unité) : Rapport d’ amortissement. 

a> 0 
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Nous obtenons donc une équation différentielle linéaire du second ordre à coefficients constants avec 
second membre. 


III.2 Solution de l’équation différentielle du mouvement 

La solution générale de cette équation différentielle est la somme de deux termes : 

• Une solution de l’équation sans second membre : solution homogène x H (t), 

• Une solution de l’équation avec second membre : solution particulière x P (t ). 

La solution totale de l’équation du mouvement sera donc : x(t) — x H (t) + x P (t) 


111.2.1 Solution homogène : 

La solution homogène correspond à la solution de l’équation différentielle sans second membre : 

x + 28 x + ù) 2 0 x — 0 

Il apparaît que la solution de l'équation différentielle homogène est tout simplement la solution trouvée 
pour l'oscillateur harmonique amorti en régime libre dans le cas des oscillations faiblement amorties : 

x H (t ) = C e~ St sin(m a t + 0) avec oi a — yjui 0 2 — S 2 

Remarque : 

• La solution générale de l’équation sans second membre 
correspond à un régime transitoire (qui ne dure qu’un 
certain temps). 

111.2.2 Solution particulière : 

Lorsque la composante x H (t) devient vraiment négligeable, il 
ne reste plus que la solution particulière, qui est la solution 
imposée par la fonction d'excitation. Nous disons que nous sommes en régime forcé ou régime 
permanent. 

La force excitatrice oblige le système mécanique à suivre une évolution temporelle équivalente à la 
sienne. Donc si F ext est une fonction sinusoïdale de pulsation m; alors la solution particulière x P (t) sera 
une fonction sinusoïdale de même pulsation co. 

Les oscillations de la masse ne sont pas forcément en phase avec la force excitatrice et présente un 
déphasage noté cp. La solution particulière correspondant au régime permanent s’écrit dont : 

x P (t)= A sin ( oit -I- (p ) 

Pour des raisons pratiques, il est commode d’utiliser la notation complexe. La grandeur complexe 
associée àjc(f) s’écrit : 

3p(t) = A e it “ t+<p) et F ext = F„e iü>t 

Déterminer les grandeurs A et tp revient à chercher le module de l’amplitude complexe. 



III.2.2.1 Calcul de l’amplitude A 

3 P (t) Vérifie l’équation différentielle avec second membre : 3‘p + 28 3 P + 6O 2 0 3p = ~ e)ut — 5e )a)t (*) 
Calculons la dérivée première puis le dérivé second : 


*(Ï) =A e j(“t+<p) ^ f SpCO = A jo> eK“ t+ <P) = jo> 3 P (t) 
3P W lsp(t) - Aj 2 o) 2 e^ a)t+<p ) = -a) 2 3p (t) 
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On remplace dans (*) et on trouve : — m 2 3 P (t) + 25jm 3 P (t) + m 2 0 3 P (t) = fle )a)t 

=^> [(m 2 0 — m 2 ) + 25mj] 3 P (t) = [(m 2 0 — m 2 ) + 25mj] A e^ < * rt+<p ) = Be ^ 

=> [(m 2 0 - m 2 ) + 25m;] Ae™ = B 

On divise sur "e- ,<p " et on trouve: [(m 2 0 — m 2 ) + 25m;] A — Be~j (p (1) 

Le conjugué de cette équation est la suivante : [(m 2 0 — m 2 ) — 25m;] A — Be w (2) 


(1) X (2) =» A 2 [(m 2 0 - m 2 ) 2 + (25m) 2 ] = B 2 => A = 

ffl.2.2.2 Calcul de tp 

[(m 2 o — m 2 ) + 25m;] A — | 


A — 

B 

< 

^ (ù) 2 q— ù) 2 ) 2 + (2S(o) 2 


= cte 


Be iv (A(m 2 0 — m 2 ) = B coscp 

<( )> < 

B (coscp — ; sintp ) { 2 ôtoA — —B sincp 


tgtp = 


-25m 


(m 2 0 - m 2 ) 


cp — Arctg (■ 


—25m 


Donc : x P (t)= 


B 


(a> 2 0 -a> 2 ) +(2 Sto) 2 


(ü) 2 0 -oo 2 ) 

mt + Arctg , 

(<w z o-*> ). 


= 5m 


-25ro 


Remarques: 

• La solution générale de l’équation différentielle s’écrit : x(t) — x H (t) + x P (t). 

• x H (t) est appelée solution homogène caractérisant un régime transitoire qui disparaît 
exponentiellement avec le temps. Quand le régime transitoire disparaît : x(t) « x P (t ) 

g 

• x P (t) est appelée solution particulière d’amplitude A — , caractérisant un 

J(lD 2 0 -lr> 2 ) 2 +(2SM) 2 

régime permanent (stationnaire) car il subsiste aussi longtemps que la force extérieure ( F ext ) est 
appliquée. Nous notons la dépendance de l'amplitude A de la pulsation m. 

• La solution x(t) aura donc souvent une allure caractéristique comme celle présentée sur la figure 
ci-dessous : 

S 


s(t) = 5 exp(-0,2t) cos(1 5t) + cos(1 6t) 



-4 ► -4 

Régime transitoire Régime permanent 
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CHAPITRE III 


Oscillations forcées amorties 


Système à un degré de liberté 2010-2011 


III.3 Etude du régime permanent : phénomène de résonance en amplitude 

III.3 .1 La variation de l’amplitude en fonction de la pulsation de la force pour différentes valeurs de £ 

Soit A(œ) l’amplitude de la solution particulière caractérisant le régime permanent (forcé) : 


d(m) = 


B 


d(m) = 


B 


7(m 2 0 - m 2 ) 2 + (2 8cS) 2 

B/ù ) 2 0 




(Jü 2 0 


N 


(l - 4 -) +( 2 «) 2 (^-) 2 ( l -(— ) 2 ) +(— ) 2 (— ) 2 ( l-(— ) 2 1 +( 2« 2 (— ) 2 

V 6) z o yj\ "0 / Vo/ "0 A \ w 0 / "0 

Fo 


s 


avec B — 


F o 


Tels que : A n — et £ = — 

é*) z o 60 0 m 


d 0 = -^V = y Donc : A 0 = Ç 

méo z o k k 



A 

lo 

1 

'(- O 2 ) 

( + (2fV(^) 2 


Posons r = — =» rl(r) = 


■ et cherchons la valeur maximale de d(r). 


l(l-r 2 y+(2fr) 2 

A(ù) ) est maximale quand le dénominateur est minimal. 

Posons r = — =» A(r ) = (1 — r 2 ) 2 + (2^r) 2 et cherchons la valeur maximale de rl(r). 
o> 0 


•d(r) max 
cLA(r) 


2^2 

cLA(r) 

dr 


= 0^ [2(1 - r) 2 (— 2r) + 8^ 2 r] = 0 


dr 


= 0 


-4r[l - r 2 - 2^] = 0 => { r J 


On a un minimum ou un maximum selon le signe du deuxième dérivée : 


d 2 A(r ) 
dr 2 


— 12r 2 + 8^ 2 — 4, donc : 


r = 0 => 


d*A(r) _ oc 2 


dr 2 


- 4 < 0 ; car 0 < £ < 1 


;r = Jl - 2^ 2 => = 8 - 16ij 2 >0; car 0 < <( < 1 

• O) 

— =0 solution refusée 

«o 


O) 


(On 


= Vl-2^ 2 

Donc : est maximale pour = ^/l — 2^ 2 < 1 ^ 

La variation de l’amplitude en fonction de la pulsation de la force pour différentes valeurs de ^ est 
représentée sur la figure suivante : 

d(co) ! 



O) 

ù> 0 


Inefficace 


Efficace 
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Remarques : 

• L’amplitude A augmente quand le rapport d’amortissement ^ diminue. 

• L’amplitude de vibration atteint un maximum quand co = cô 0 : on dit qu’il ya résonance : la 

valeur maximale ( — J (correspondant à la valeur maximale de l’amplitude) n’est pas égale à 1 
\ù) 0 y 

mais égale à yj 1 — 2Î, 2 < 1 . 

Discussions: 

> Au début d’un mouvement résonnant, lorsque la force est appliquée au système, la majeure partie de 
l’énergie, fournie lors de chaque cycle, est emmagasinée dans le système; une faible partie se dissipe en 
frottement. L’énergie ainsi emmagasinée par le système fait augmenter progressivement l’amplitude de ses 
oscillations jusqu'à une valeur maximum. Cette valeur subsiste tant que subsiste l’apport d’énergie par la 
force extérieure. 

> Plus l’amortissement est faible, plus cette courbe est aigue et plus le maximum est grand : en l’absence 
d’un amortissement suffisant rien ne viendrait limiter les amplitudes des oscillations à s’amplifier, risque 
de destruction du système : le système entre en résonance. Les conséquences peuvent être graves. On peut 
citer deux cas connus : 

CE) Le 18 avril 1850 à Angers, un régiment traversant au pas cadencé (harmonieux) un pont 
suspendu enjambant le Maine provoqua sa destruction. 

CE) Le 7 novembre 1940, six mois après son inauguration, le pont suspendu de Tacoma (Etats- 
Unis) était détruit par les effets des rafales de vent qui sans être particulièrement violentes 
(60 km.tf ! ) étaient régulières. 


III.3 .2 La variation de la phase en fonction de la pulsation de la force pour différentes valeurs de £ : 

Soit (p la phase initiale de la solution particulière caractérisant le régime permanent (forcé) tel que : 


tgtp = 


-2ôo) 

(' ü) 2 0 -(O 2 ) 


• Nous remarquons que tg <p est négatif. Cela paraît normal qu'il y ait un retard de l'oscillateur par 
rapport à la force qui entretient le mouvement. L'oscillateur harmonique essaie de suivre le 
mouvement en étant ralenti par les frottements, donc il doit obligatoirement prendre du retard par 
rapport à l'oscillation excitatrice donc avoir un déphasage négatif. Ce déphasage est dépendant de 
la pulsation de la force ai. 


-2ÔO) -2ÔCJ 



• L’oscillateur est toujours en retard de phase par rapport à la force et ce retard augmente lorsque la 
pulsation augmente. 
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• L’oscillateur et la force sont en phase pour O) — 0 

Conclusions : 

Selon la valeur de ^ on a 3 cas possibles : 


1 er cas : Faibles fréquences : ^ « 1 (o> « o> 0 ) -> 
2éme cas : Hautes fréquences : ^ » 1 (o> » o> 0 

3éme cas : La résonance : ^ = 1 (o> « o) 0 ~ o> r ) -> 


A ~ An = T 


(p — 0 

{ a »° 
W = - 


n 


(-)* = vi-2ç 2 , 

\ù)oJ Aq 


n 

V = -2 


Remarques: 


î 

2C/1-S 2 . 


• Si ^ = 0 ( système non amorti) : l’amplitude tend vers l’infini or en réalité, les systèmes sont tous 
amortis donc l’amplitude n’est jamais infini. 

A 1 

• Si ^ « ( système faiblement amorti) : - n} ^- ~ ; (o> « o> 0 « o> r ) 

III.3 .3 Phénomène de résonance et Facteur de qualité 


• Dans les systèmes éléctriques, ce phénomène permet de calculer le facteur de qualité Q qui 

Amax ^ t 
A o ~ 2f‘ 

üio 


augmente lorsque l’amplitude maximale augmente Q — 


Une autre méthode pratique pour déterminer le facteur de qualité : Q 


ù>2-(Oi 


• Pour caractériser l’acuité (intensité) de la réponse d’un oscillateur en fonction de la pulsation, on 
définit une bande passante : (o 2 - (o i 



Conclusions : 


• Quand ï, augmente => Q diminue ^ o ) 2 — oj 1 augmente => la courbe de résonance est plus 
large => diminution de l’amplitude de résonance donc de la qualité aussi. 

• Les extrémités de la bande passante correspondent à une amplitude de vitesse sfî fois plus petite 
qu’à la résonance. 
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B3 Points clefs 
Oscillations forcées amorties 


1. Un mouvement unidimensionnel x : 


L’équation de Lagrange : 

Jt ê) 'Tx = ~ S +F ext ; D = 7 et x 2 , F ext = F 0 sin ait 


dx 


l'équation différentielle sous la forme réduite : x + 28 x + a) 2 0 x = ^sin ait 
Solution de l’équation différentielle du mouvement: x(t) = x H (t) + x P (t) 

La solution homogène: x H (t ) = C e~ St sin(m a t + 0) avec iv a = ^/w 0 2 — S 2 
La solution particulière : x P (t)= , === , 2 sin (ait + Arctg ( ~ 2g “ ) 

La variation de l’amplitude en fonction de la pulsation de la force : 




( 1_ Q 2 ) +(2f)2 Q 2 


max ^ ma ^ e pour = ^/l — 2^ 2 




= 1 


= - c 




Ç = 0 => = 0 <=>(£> = 0 ou <p = — 7T 

Fo 


> Faibles fréquences : ^ « 1 (o> « o) 0 ) 

> 


Ao = f 


<p = 0 

t A ~ 0 

Hautes fréquences : ^ » 1 (<o » (o 0 ) | _ _ jr 

( — li _ 7?2 _ 1 

w v 1 ^ * a 0 2 fjï=^ 


> La résonance : ^ = 1 (<o » o) 0 ~ to r ) 


n 

( P = ~ 2 


Le facteur de qualité : Q = 


2. Un mouvement rotationnel 0 : 


L’équation de Lagrange : 


ù) 0 


2f (»>2 (»)] 


rf 

dt 


(S) l Fe * f l ’ r : est la direction d’action de la force F ext 
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CHAPITRE IV 

Oscillations libres des systèmes à plusieurs degrés de liberté 

Introduction : Dans ce chapitre, nous examinons les systèmes qui se composent de deux ou plusieurs 
oscillateurs qui sont couplés dans une certaine façon et qui ont plus d'une pulsation d'oscillation. Nous 
allons voir que ce couplage produit de nouveaux et d’importants effets physiques. Chacune des pulsations 
correspondent à une manière différente dans laquelle le système peut osciller. Ces différentes façons sont 
appelés « modes normaux ». Les modes normaux d'un système sont caractérisés par le fait que toutes les 
parties du système oscillent avec la même pulsation. Les oscillateurs sont couplé parce qu’ 
ils se trouvent rarement dans un isolement complet et sont généralement capables d’osciller avec de 
différentes façons. Les oscillateurs couplés sont également importants car ils ouvrent la voie à la 
compréhension des ondes dans les milieux continus. Le mouvement des ondes dépend des systèmes 
voisins qui vibre et qui sont couplées entre elles et peuvent donc transmettre de l’énergie entre elles. 
Définition : 

Un système est à plusieurs degrés de liberté (ddl) si plusieurs coordonnées indépendantes sont 
nécessaires pour décrire son mouvement. Il y a autant d’équations de Lagrange que de degrés de liberté 
ou de coordonnées généralisées. 

IV. 1 Systèmes à 2 degrés de liberté 

Pour l’étude des systèmes à deux degrés de liberté, il est nécessaire d’écrire deux équations différentielles 
du mouvement que l’on peut obtenir à partir des équations de Lagrange : 

{ d / dL \ / dL \ 

dtvdçq/ KdqJ 
d / dL \ / dL \ 

dt \dq 2 ) \dq 2 ) 


Un système à 2 degrés de liberté possède 02 coordonnées généralisées, 02 équations différentielles et 02 
pulsations propres ( m 1; o> 2 ). 


IV. 1 .1 Les type de couplages 


a) Couplage Elastique : Le couplage dans les systèmes mécaniques est assuré par élasticité. Dans les 
systèmes électriques, on trouve les circuits couplés par capacité, ce qui est équivalent au couplage 
par élasticité. 



Les équations différentielles correspondantes sont : 


x 1 + 2ô 1 x 1 + eo 2 x 1 = a ± x 2 

, Tels que : a 1 x 2 et a 2 x ^ sont les termes de couplage. a i et a 2 sont des 

x 2 + 2 S 2 x 2 + O) x 2 = a 2 x 1 

constantes. 
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b) Couplage Visqueux : Le couplage dans les systèmes mécaniques est assuré par amortisseur. Dans 
les systèmes électriques, on trouve les circuits couplés par résistance, équivalents au couplage par 
amortisseur. 



Les équations différentielles correspondantes sont : 


pîq + 2S 1 Âq + to 1 2 3 4 5 x 1 = b 1 x 2 
uiq + 2 S 2 x 2 + ù) 2 x 2 = b 2 x i 


Tels que : b 1 x 2 et b 2 x 1 sont les termes de couplage. b 1 et b 2 sont des 
constantes. 


c) Couplage Inertiel : Le couplage dans les systèmes mécaniques est assuré par inertie. Dans les 
systèmes électriques, on trouve les circuits couplés par inductance, équivalents au couplage par 
inertie. 




Les équations différentielles correspondantes sont : 

Çx 1 + 2S 1 x 1 + o) 2 x 1 = c 1 x 2 

, , Tels que : c 1 x 2 et c 2 x 1 sont les termes de couplage. c 1 et c 2 sont des constantes. 

Lx 2 + 2S 2 x 2 + ù) x 2 = c 2 x 1 


IV. 1 .2 Méthode générale de résolution des équations de mouvement. 

Pour un système mécanique, la mise en équation du système couplé passe par la méthode à suivre 
suivante : 


1 - On écrit les 2 équations différentielles en fonction des coordonnées généralisées. 

2 - On fait l’hypothèse que le système admet des solutions harmoniques. Ce qui signifie que le système 
peut osciller avec la même pulsation pour tous les oscillateurs. 

3 - La résolution des systèmes d’équations permet d’obtenir 2 pulsations particulières œ 1 et œ 2 ; 

ce sont les pulsations propres. 

4 - On substitue ensuite œ 1 dans l'une des 2 équations et l’on obtient le I e1 mode propre. On substitue 
ensuite œ 2 dans l'une des 2 équations et l’on obtient le 2 eme mode propre. 

5 - On écrit les 2 solutions générales des équations différentielles du mouvement. 
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IV.l .3 Exemples de systèmes a 2 DDL 

IV. 1 .3.1 Pendules couplés : (Couplage Elastique) 

Considérons deux pendules qui sont couplés par un ressort horizontal de 
constante de raideur k à une distance a de l'axe de rotation. 

£0 1. Equations différentielles du mouvement : 

❖ Les coordonnées des éléments du système : 

La masse m 1 se trouve à une distance de O. 


x. 


nij 


m-y 


— l v sin 9 ± 


x. 


m ! 


= Z^cos 9 ± 


[ym 1 - ~k- cos 9 1 [ÿ mi = l 1 0 1 sin 9 x 
La masse m 2 se trouve à une distance Z 2 de O. 


17 2 _ / 2p2 

v ml ~ L 1 



m, 


m 2 


x m7 = l 2 .sin 9 2 


x m? = co s 9 2 


[ym 2 - -h- cos 9 2 [ÿ = l 2 è 2 sin 9 2 


v‘ 


m 2 


= h 2 '0\ 


k — { a.sinOx — a.sin0 2 — a(sin0 1 — sin0 2 ) 

❖ L’énergie cinétique du système : T = Tm 1 + Tm 2 = \ mv m 2 + \mv m 2 


Tm 1 = -m 1 (x rni +ÿ m J 


2 — (cos 9 X 2 + sin 9 1 2 ) 


Tm 1 = \m x l 2 x G\ 


Tm 2 = -m 2 (x m2 + ÿ m2 ) 2 = ^m 2 l 2 2 9i(cos 9 2 2 + sin 9 2 2 ) 


Tm 2 = -m 2 l 2 2 G 2 


=> T = - ( m 1 l 2 1 0l + m 2 l 2 2 G 2 2 ) 

❖ L’énergie potentielle du système : U — U k + u mi +u m2 

Si on choisi comme origine des énergies potentielles l’axe (Ox) on a pour les deux masses : 

U mi +U m2 — cos G 1 —m 2 gl 2 . cos G 2 (Le signe moins vient du fait que la masse m est inférieur 

à l’axe choisi). 

1 

U = -ka 2 (sinG l — sin6 2 ) 2 —m 1 gl 1 . cos G l —m 2 gl 2 . cos G 2 


■ La fonction de Lagrange sera donc : 

L = T-U = \m x l 2 j\ + \ml 2 2 Q\ - \ 


ka 2 (sinG 1 — sinG 2 ) 2 +m 1 gl 1 . cos G 1 +m 2 gl 2 . cos G 2 


On remarque bien deux coordonnées généralisées qui décrit le mouvement donc on aura deux équations 
de Lagrange : 


d ( dL ^ ( dL} n 

dtvaej VaaJ - 0 

d ( dL ^ ( dL\ n 

dt\dê 2 ) \d0 2 )~ ° 


^ — ( -—7-\ — 


dt \d9 ly 
/ dL \ 

—— 1 = —ka^ cos0 1 (sin0 1 — sinG 2 ) — m 1 gl 1 sinG 1 

\Ou -jy 

=> m 1 l 2 1 G 1 + ka 2 cosG 1 (sinG 1 — sinG 2 ) + m 1 gl 1 sinG 1 = 0 

( sin 9 ~ 9 


Dans le cas des faibles oscillations, les angles sont très petits on a : 


cos 9 


i-£ 

2 
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=> m 1 l 2 1 G 1 + ka 2 (G t - G 2 ) + = 0 

Donc les 02 équations différentielles du mouvement sont : 

(m 1 l 2 1 Ô 1 + ( ka 2 + m 1 gl 1 )G 1 = ka 2 G 2 (1) 

[m 2 l 2 2 G 2 + (ka 2 + m 2 gl 2 )G 2 = ka 2 G l (2) 

Remarque 

• Le terme de couplage ka 2 est en fonction de k donc le couplage est élastique. 

• Si a — 0 ou k = 0 => couplage nul : les deux systèmes sont indépendant. 

• Les deux équations différentielles possèdent 02 solutions 0 1 (t)et G 2 (t). 

£Q 2. On fait l’hypothèse que le système admet des solutions harmoniques : 

Donc : G ± (t) — A^inÇcot + (p) et 0 2 (t ) = A 2 sin(cot + tp 1 ) 

Tels que : A lf A 2 , tp et tp', to est l’une des pulsations propres du système. 

( G ± (t) — A 1 sin(ù)t + tp) => G 1 — —oi z Q 1 
{0 2 (t) = A 2 sin(cot + tp 1 ) =3 Ô 2 = —(jô 2 0 2 
On remplace dans les équations (1) et (2) donc : 

((ka 2 + m 1 gl 1 - m 1 l 2 1 to 2 )G 1 - ka z G 2 = 0 (3) 

{ -ka z 6 1 + (ka 2 + m 2 gl 2 - m 2 l z 2 (o 2 )G 2 = 0 (4) 

£Q 3. Calcul des pulsations propres : On suppose que m! = m 2 = ni, l x — l 2 — l 
(ka 2 + mgl - ml 2 a ) 2 — ka 2 \(o 1 \_ro\ 

\ —ka 2 ka 2 + mgl - ml 2 co 2 ) 


Ces deux équations accepteront une solution si le déterminant =0 

ka 2 + mgl — ml 2 œ 2 —ka 2 

—ka 2 ka 2 + mgl - ml 2 co 2 


= 0 (ka 2 + mgl — ml 2 co 2 ) 2 — (ka 2 ) 2 = 0 


2 2 

(ka 2 + mgl — ml 2 co 2 ) - (ka 2 ) = 0 => ka 2 + mgl — ml 2 


(O 




/» 2 - 9 
"2-y 


tels que : j 


o) 1 : la 1ère pulsation propre 
to 2 : la 2ème pulsation propre 


+ka 2 
— ka 2 


Remarque 

• Si a — 0 ou = 0 , le couplage est nul => (n 2 1 = (o z 2 — j 

• Lorsque le système oscille avec une de ses 02 pulsations on dit que le système oscille dans un de 
ses deux modes. 

£0 4. Les modes d’oscillations 


Le mode c’est l’état dans lequel les éléments dynamiques du système effectuent une oscillation 
harmonique avec la même pulsation qui correspond à une de ses deux pulsations. 

4.1 Calcul des modes d ’ oscillations : 


Dans chaque mode les deux masses effectuent des mouvements harmoniques simples avec la même 
pulsation (<*>! ou to 2 ) et les deux pendules passent par la position d’équilibre au même instant. 

Premier mode : on remplace dans (3) ou (4) par co 2 1 — j + 2 (y) 2 : 

On obtient après calcul : G 2 = —G { 

Remarque : 

• Dans le premier mode les deux pendules ont la même pulsation oq , la même amplitude et un 
déphasage n . 
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• Les deux pendules ont des mouvements opposés. 

• Elongation et compression du ressort chaque période sauf au point du milieu du ressort. 



Deuxième mode : 

on remplace dans (3) ou (4) par o) 2 2 — j ■ 

On obtient après calcul : G 2 = 0 x 

Remarque : 

• Les deux pendules se déplacent dans le même sens. 

• Le ressort ne subit aucune variation de sa longueur. 




CQ 5. Calcul des solutions des équations différentielles : 

Chacune des mouvements 9i et 02 possède deux composantes harmoniques de pulsations oq ou oq 
Comme les équations différentielles sont linéaires, toute combinaison de solutions reste solution du 
système. 

La solution générale s’écrit alors comme une combinaison linéaire des deux solutions. 

— A^inÇütit + q){]+B 1 sin{o) 2 t + <p 2 ) 

1 0 2 (t) — i4 2 sin(m 1 t + tp{)+B 2 sin{o) 2 t + tp 2 ) 

Dans le premier mode : co = o) 1 => 0 2 = ~8\ =* A± = —A 2 =» L 1 (_ 1 1 ), l^est le 1 er vecteur propre 
Dans le deuxième mode : co = to 2 =^9 2 = 6 1 ^B 1 — B 2 ^> F 2 (*), L 2 est le 2 eme vecteur propre 
Donc : 

( 0\(t) — AsinÇoiit + tp 1 ) + Bsin((o 2 t + tp 2 ) 

(0 2 (t ) = — AsinÇo)^ + tp 1 ) + Bsin(ü) 2 t + tp 2 ) 


6. Calcul des constantes A. B. tp 1 et tp 2 


Supposons que 


|0 1 (t) = Oo.Ô 1 (t) = 0 

[.0 2 (t) = O,Ô 2 (t) = O 
\ 0l(t) = AcdiCOS^Wit + + B(j0 2 C0S(C0 2 t + (p 2 ) 

[9 2 (t) — —AùiiCOsÇcôît + (p t ) + Bù) 2 cos(ù) 2 t + tp 2 ) 


C 0i(O) = Asin(jtp{) + Bsin(tp 2 ) = 6 0 
(0 2 (O) = —Asin(<p{) + Bsin(tp 2 ) = 0 


et 


0i(O) = Aw^ositpi) + Ba) 2 cos((p 2 ) 
[Ô 2 (0) = -AMiCOsitpi) + Bù) 2 cos(tp 2 ) 
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Donc : 


n 


<Pi = <P2 = ±2 



9i (t) = y [ sin (wjt + 0 + sin{o) 2 t + ^)] 
0 2 (0 = y [-sin («it + f) + sin(o) 2 t + 1)] 


01 (t) = OqCOS ( ^ 2 - — ) t + COS ( — y Ml ) t 

0 2 (t) - -0 o sin ( — - — ) t + sin ( — y ^ t 
CQ 7. Phénomène de battement : 

Lorsque le couplage est faible (k faible), les pulsations propres des 2 oscillateurs (uq et o> 2 ) sont voisines 
(<*>! « <» 2 =» Am = <n 2 — U) 1 est faible), il se produit un phénomène de battement. Les 2 oscillateurs 
se transmettent de l’énergie entre eux et vibres avec une pulsation œ égal à la moyenne des deux 

pulsations propres a) — - (a) 2 + cof) avec une période égale à T — — — ^ +ù) Tandis que la pulsation 

1 4-7T 

du battement est égale à a) B — - (<w 2 — (of ) , avec une période T B — — — — 



IV. 1 .3.2 Pendules couplés : (Couplage inertiel) 

Considérons deux pendules qui sont couplés par une masse m 1 qui 
se trouve à une distance f de l'axe de rotation. 


1. Equations différentielles du mouvement : 

❖ Les coordonnées des éléments du système : 

La masse m 1 se trouve à une distance Z x de O. 

\ x m± = l^sin 6 1 (x mi = l^cos 9 ! 

I ym x = ~k- cos 6 i [ÿ mi = Z À sin G ! 

La masse m 2 se trouve à une distance (f + Z 2 ) de O. 


ni] 


v 


ml 



= i/e 


x. 


m 2 


Calcul de v 


m 2 


— Z x . sin 6 1 + l 2 .sin 8 2 


Xm 2 — lidiCOS 8 1 + Z 2 0 1 cos 8 2 


[ym 2 — ~h- cos 0 i — l 2 - cos 0 2 [ÿm 7 — W&i sin 8 x + l 2 d 2 sin 8 2 


m 2 


■v 1 2 m 2 — (Zi0icos 8 ! + l 2 e 2 cos 8 2 ) 2 + (l 1 8 1 sin 8 1 + l 2 d 2 sin 8 2 )‘ 


v z m — l z 1 8 2 + l z 2 d 2 + 2 Z 1 0 1 cos 8 l 2 d 2 cos 8 2 + 2 1 1 8 1 sin 8 x . l 2 d 2 sin 8 2 


û2 
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v 2 m2 — Z 2 X 0 2 + Z 2 2 0f + 2Z 1 0 1 Z 2 0 2 (cos 0 i cos 0 2 + sin 0 x sin 0 2 ) 

Or : cos (0 i — 0 2 ) = cos 0 -^os 0 2 + sin0 1 sin 0 2 ) => v 2 mz = Z 2 i0 2 + Z 2 2 0f + 2 Z 1 0 1 Z 2 0 2 cos (0 x — 0 2 ) 

Dans le cas des faibles oscillations, les angles sont très petits on a : {cos (0 x — 0 2 ) » 1 
Donc : v m2 — l ^0 2 + Z 2 0f + 2Z 1 Z 2 0 1 0 2 = + Z 2 0 2 ) 2 

❖ L’énergie cinétique du système : T = Ttrix + Tm 2 = ~ L mv mi 2 + ^mv m2 2 
' Tm 1 = - 2 m x {x mx + ÿmj 2 = ^m 1 Z 1 2 0 1 2 (cos0 1 2 + sin0 1 2 ) 

=> Tm 1 = 

■ Tm 2 = \m 2 v 2 m2 = \m 2 (l 1 é 1 + Z 2 0 2 ) 2 

7 = in^Z 2 ^ + im 2 (Z 1 0 1 + Z 2 0 2 ) 2 

❖ L’énergie potentielle du système : U — U mi +U m2 

Si on choisi comme origine des énergies potentielles l’axe ( Ox ) on a pour les deux masses : 

U mi +U m2 — —Hiygly. cos 0 1 —m 2 g(l 1 . cos 0 1 + Z 2 . c os 0 2 ) 

(Le signe moins vient du fait que la masse m est inférieur à l’axe choisi). 

U = —gl\(m 1 + m 2 ) cos G 1 —m 2 gl 2 . cos 0 2 

■ La fonction de Lagrange sera donc : 

L = T — U = — mil 1 0 1 + — m 2 (Zi0i + l 2 0 2 ) + gliipii + m 2 ) cos 8 l +m 2 gl 2 . cos 0 2 


On remarque bien deux coordonnées généralisées qui décrit le mouvement donc on aura deux équations 
de Lagrange : 



— + m 2 Zi(Z 1 0 1 + Z 2 0 2 ) 

=> m 1 l 2 1 G 1 + m 2 l 2 1 G 1 + m 2 l 1 l 2 G 2 + gli(m 1 + m 2 ) G 1 = 0 

-) = ~gl i("k + m 2 ) sin 0 X 

L 2 

=> (m 1 +m 2 )Z 2 0i + gl 1 (m 1 + m 2 ) G t = - m 2 l 1 l 2 G 2 
On divise sur (mi+mz^i et on trouve : l 1 G 1 + gG t = — - — mi 2 0 2 

(mi + 1112)11 



fd/5L\ 

dïfej _m2f2(fl01 + f202) 

( dL\ 

{ w J = -^oh-^e 2 


Ul 2 Z 2 (^1®1 + Z 2 0 2 ) + m 2gW- 8 2 — 0 


=> m 2 l 2 G 2 + m 2 gl 2 . G 2 — —m 2 l 2 l^G ^ 


On divise sur m 2 Z 2 et on trouve : Z 2 0 2 + tj0 2 = — Z]^ 


IiÔi + *®i = -^^Ô 2 (1) 

^2®2 + 002 — — ^1^1 (2) 

£Q 2. On fait l’hypothèse que le système admet des solutions harmoniques : 

Donc : 0i(t) = A^inÇcot + <p) et 0 2 (t) = A 2 sin((ot + tp 1 ) 


Donc les 02 équations différentielles du mouvement sont : 
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Tels que : A t , A 2 , tp et tp', oi est l’une des pulsations propres du système. 

( 0i(t) = A^inÇoit + tp) =$ = —oi 2 0 1 

{e 2 (t) — A 2 sin{oit + tp 1 ) =» 0 2 — —oi 2 G 2 

On remplace dans les équations (1) et (2) donc : 

o ( 3) 

(O - i 2 m 2 )g 2 - i 1 a) 2 e 1 = o (4) 

CQ 3. Calcul des pulsations propres : On suppose que m 1 - m 2 — m,l 1 — l 2 — l 

{ {g-loi 2 )0 1 - l -oi 2 0 2 = 0 (4) 

((g — loi 2 )G 2 — loi 2 G — 0 (5) 


(g -loi 2 

loi 2 \ 

~ K?) 

(N 

3 

1 

g - loi 2 ) V 2 


Ces deux équations accepteront une solution si le déterminant =0 

loi 2 


g-UA - — 
-lot 2 g -loi 2 


f 2 \ 2 1 s ^ 

= 0 [g — la) ) — - ( la) ) = 0 ; C’est l’équation aux valeurs propres. 


Oi 2 1 


« 1+ vl> 


Oi 2 2 = 


{g -loi 2 ) --{loi 2 ) =0 =$g-loi 2 = 


Ce sont les valeurs propres. 


r 1 7 

+ — loi 

V2 

1 2 

—loi 2 

V2 


oii valeurs propres. 




3.Calcul des modes d’oscillations ou les vecteurs propres : 


Premier mode : on remplace dans (5) ou (6) par oi 2 1 — — — j- 

^ 1+ V2-* 

K S -t~) 

„ K1 + -S) 

(9 - K — = 0 

<(1 + ^ 


(1- j-p-) G t — t - 0 2 => 0 2 — V2 G x : c’est le premier mode 


V2 


2 ^ 1+ vf ) 


g 


Deuxième mode : on remplace dans (5) ou (6) par oi 2 — ± 

V5-* 

On trouve : 0 2 — — \[2 O^. c’est le deuxième mode. 

£Q 4. Calcul des solutions des équations différentielles : 


Chacune des mouvements 9i et 6i_ possède deux composantes harmoniques de pulsations uq ou to 2 . 
Comme les équations différentielles sont linéaires, toute combinaison de solutions reste solution du 
système. 

La solution générale s’écrit alors comme une combinaison linéaire des deux solutions. 

f (t) = AV^sinÇoi^ + + E?V2sin{oi 2 t + tp 2 ) 

l0 2 (t) — AVïsinioiît + tp 1 ) — BV2sin(oi 2 t + tp 2 ) 
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IV.l .4 Généralisation aux systèmes à n degrés de liberté : Principe des opérateurs 
IV.l .4.1 Energie cinétique généralisée (Opérateur associé à l’énergie cinétique ) 

Un système à n degrés de liberté possède n variables : x lt x 2 , x 3 , x n . 

*] O 1 O 1 O 

L’énergie cinétique généralisée ; T = -m 1 x{ + -m 2 x 2 + ••• - m n x„ 

Sim! — m 2 — ■■•m n = m=> T = ^m(xf + xf + ••• x =» T = îxf (1) 

/. 

Soit le vecteur vertical (colonne) : |x) = 

\-^/i 

Soit le vecteur horizontal (ligne) :(x| = (x 1( x 2 , ... x n ) BRAS 

l.\ L 


X 1 I 

*2 KIT 


<*!*> = (x ll x 2 ,...x n ) 


Xl 


T = -(x 1 ,x 2 ,x 3 , 


x ? ] = £f=l Xi <=> (x\x) = ( X lt X 2 , ... X n ) 

\ x n 


• x n ) 


xi 


*2 = Zi= ixf => T = -m(x\x) 


[m 

0 

0 

■■ 0 \ 

/x i\ 

0 

m 

0 

0 


f x 2 \ 

0 

0 

m ■■■ 

0 


* 




0 


J / 

\o 

0 

0 ... 

... m) 

\ n / 



T : Matrice carrée (nxn), Opérateur associé à T. Les éléments de T sont déduits des dérivées 
IV.l .4.2 Energie potentielle généralisée (Opérateur associé à l’énergie potentielle) 

T=;Mi + + -\Kxi • 


Si A 


U = \k{x \x) =» U = \{x\U\x ) . Tel que : U = 


Xn) 

=> 

U = 

i*j 

f k 

0 

0 

°\ 

0 

k 

0 

0 

0 

0 

k 

0 




0 

\o 

0 

0 

fc/ 


( 2 ) 


T l : Matrice carrée (nxn), Opérateur associé à U. Les éléments de TL sont déduits des dérivées ^ 

IV.l .4.3 Equation différentielle 

Le Lagrangien (la fonction de Lagrange) : L — T — U — ^(x\T\x) — ^(x\ZL\x) 

Equation de Lagrange (£) = 0 

dL 1 â 

ox t 2 dx t 

.dL 1 d 

k— = -~-r-{x\U\x) 

vdxi 2 dx t 

On a : -^-(x|T|x) = 2(Ii\T\x) et -^-(x|Xt|x) = 2{Ii\TL\x) avec I t vecteur unité. 

o Xi u x î 

(h\ = (1, 0,0,0 0), </ 2 1 = (0,1 AO 0) (I n I = (0,0 ,0,0 1). 

dL d (dL\ 

— = (/,m*>^-(— )= 


Oï 

Oj 

O) 

O 

O. 
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= —{Ii \V.\x) => (Ii\T\x) + (Ii\V.\x) — 0, c’est l’équation de Lagrange. 

~ rm + «l*> = » - l*> + r‘W = o =» ^ w + xw = o 

£ — T 1 .U : Opérateur associé au Lagrangien L, T -1 : Matrice inverse de T . 


IV.l .4.4 Equation aux valeurs propres : 

d 2 

La solution de l’équation —, |x) + L\x) — 0 peut être sous la forme complexe : Ce J0)t 


- — (*) 2 \x) => — (*) 2 \x) + £\x) = 0 =>(£ — ù) 2 )\x ) — 0 , 


(£ — ru 2 ) |x) = 0 => Det [£ — ru 2 /] = 0, c'est l’équation aux valeurs propres, ru t : valeurs propres. 

♦> valeurs propres : 

Exemple : £ = , 

\£21 l 22/ 

T , ✓ .. , r 2i (L 11 ^12\ 2 fl 0\ / 1*11 — (A 2 £i 2 \ 

L équation aux valeurs propres : £ — <*r/= r r — ru ' I J = „ 

4 F F U 21 1*22/ Vo 1/ \ L 21 L 22 -ù> 2 J 

— 0 => (£11 - ü) 2 )(L 22 — O) 2 ) — Li2- /-2i — 0 


/-Il — < W /-12 

/-21 /-22 — É 

A chaque valeur propre tUj correspond un vecteur propre If , Wx ^ If 

Lu ru 1 2 L 


❖ Pour ; a) = o) 1 : (£ — ni! 2 ) V ± — 0 


'21 


Z-22 — 


((Lu — rUi 2 )^!! + £ 12^12 — 0 ? 

(1*21*11 - (L 22 - Wl 2 )^12 — 0 11 12 ' 


❖ / J orrr ; tu = tu 2 : (£ — to 2 2 ) V 2 = 0 


/-Il — 

L 2 i 


L 22 — w 2^ 


îte:)^~^te:) 

te:) = o 

te:) = (S) 


=> *21 et * 22 ? 

IV.l .4.5 Solution des équations différentielle : 

^iWV A îr e 7(ûJit+<pi) , By* e i(.o>2t+(P2) = A e Ko)it+(Pi ) + g e Ku 2 t+(P2) 

V*2 (t)l 1 2 l* 12 ) \ x 22/ 


(Xi(t) = Ax u e Kû>lt+ ^ + Bx 21 e^ t+ ^ 
U 2 (t) = A* 12 e Kwit+<Pl) + Bx 22 eJ^ t+( P^ 
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